Lower bound for ground state energy of BEC in a rotating optical lattice 
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We use the frustrated XY model approximation of BEC in a rotating optical lattice and formulate 
the problem of the ground state in terms of eigenvectors and eigenvalues of frustrated adjacency 
matrix (coupling matrix) . By using this formulation, we show that there is a lower bound for ground 
state energy in terms of maximum eigenvalue of this matrix. 

PACS numbers: 



I. INTRODUCTION 

Effect of the rotation of optical lattice on physical 
properties of BEC is an interesting problem. This ef- 
fect can be observed as vortex in lattice and in this way 
is related to many problems in condensed matter physics, 
e.g. Josephson junctions arrays, ... X 

It can be shown that with appropriate 
approximations^—, the problem of the ground state 
of BEC in a rotating optical lattice can be formulated 
as frustrated XY model. In this formulation, rotation 
plays the role of frustration parameter— ~—. 

Here, we focus on ground state of BEC in a rotating 
optical lattice in this approximation and by use of prop- 
erties of frustrated XY model, we find a lower bound for 
ground state energy for different optical lattices. 



II. BEC IN ROTATING OPTICAL LATTICE: 
FRUSTRATED XY MODEL 

It can be shown that in appropriate approximation, 
Hamiltoian for BEC in a rotating optical lattice can be 
reduce to frustrated XY Hamiltonian^r— , 



H = £ Ji, 



-Oj + Aij) 



(2.1) 



<i,j> 



Suppose that rii is number of atoms in ith site and Oi 
is the phase of localized wave function on ith site. Jij is 
proportional to ^fn^n] and Ay is line integral of A = VLzx 
r. If we define frustration parameter / as / = 2£lm/h, 
where m is mass of condensate and h is plank constant, 
then we can define 2irBijf — Ay for nonzero / and 



H 



J l3 cos{6 l - 6, + 27rBtf/) 

<i,j> 

can be written as 



(2.2) 



Hamiltonian [ 

H = -Re{ J2 Jijexp{i(ei - 6 3 + torB tj f))} (2.3) 



<t,3> 



where Re{x} is real part. By introducing frustrated ad- 
jacency matrix, Sij = Jijexp(i2TrBij /), we have, 



where N is number of sites and ^ = exp(—i9i)/vN. The 
factor 2 enters because we consider both and (j, i) 
terms in this case. 

Now we focus on formulating of the problem in terms 
of eigenvalues and eigenvectors of S. Suppose that r/ fc is 
eigenvector correspond to eigenvalue X k of matrix S, i.e. 
Sr/ k = X k r] k . Also we suppose that set of eigenvectors 
of S are complete and orthogonal. Then we can write 
following expansion for £, 



(2.5) 



Using this expansion in Hamiltonian 12.4] 

H = -N/2Re{^2 c k c k *X k } (2.6) 

k 

But there is some restriction on c coming from properties 
of £. Because of = 1, we have J2 k c k c k * = 1 and 
therefore \c k \ < 1 for any k. Also from definition of £j we 
have, 



which leads to, 



E 

m,k 



m k* m k* _ J_ 

c c m m — 



(2.7) 



(2.8) 



The problem of minimization of Hamiltonian 12.41 is 
equivalent to minimization of Hamiltonian 12.61 with re- 
spect to complex numbers c,-'s and with constrains! 



III. 



LOWER BOUND FOR GROUND STATE 
ENERGY 



We focus on minimization problem defined by equa- 
tions 12.61 and restrictions 12.81 

We consider simple situation, a plaquette with unit 
length and uniform Jjj, which can be shown that frus- 
trated XY model for it can be solved exactly. In 
this case, S is given as, Si 2 = 1,S23 = 1, S34 = 



exp(-2mf) 1 Sn = 1, 



and Sji — 



S*j (which is prop- 



H = -Y/2i?e{£*S£} 



erty of S). Numerical results shows that in this case, 



(2.4) E/N = -X max /2. 
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From the general Hamiltonian 12.61 and the constrain's 
on c"s (\c l \ < 1), we can see that, 

H/N > -X max /2 (3.1) 

Therefore the minimum of Hamiltonian also satisfies this 
inequality. It means that 

E/N > -X max /2 (3.2) 



IV. CONCLUSION 

We have shown that ground state energy of a BEC in 
an rotating optical lattice in frustrated XY approxima- 
tion, has a lower bound which determined by coupling 
matrix of system. Formulation of the problem in terms 
of properties of coupling matrix, maybe can be used for 
more analytical results especially relation between XY 
models with different coupling matrix. 
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